Frustrated spin model as a hard-sphere liquid 
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We show that one- dimensional topological objects (kinks) are natural degrees of freedom for an 
antiferromagnetic Ising model on a triangular lattice. Its ground states and the coexistence of spin 
ordering with an extensive zero-temperature entropy can be easily understood in terms of kinks 
forming a hard-sphere liquid. Using this picture we explain effects of quantum spin dynamics on 
that frustrated model, which we also study numerically. 
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PACS numbers: 05.50.+q, 75.10.Jm, 75.30.Fv, 64.60. Cn 

Geometrically frustrated materials recently emerged as [ 
a new broad class of solids with interesting and rather 
unusual properties |jj . While some of these systems stay 
disordered at all temperatures, others order often in an 
unexpected way, showing no universality typical for crit- 
ical behavior of conventional systems. One may won- 
der, however, to which extent the complexity of ground 
states and excitations of frustrated models is their gen- 
uine property and to which it is, essentially, a conven- 
tional behavior obscured by an unfortunate choice of vari- 
ables, in which these models are formulated. In this Let- 
ter we consider a frustrated Ising model showing a spin- 
density-wave (SDW) ordering, more common for systems 
with continuous symmetries. We explain the origin of 
this strange behavior and give a simple description of 
ground states of that model (the number of which grows 
exponentially with the system size) by mapping it on a 
hard-sphere liquid. We use this approach to study the 
role of quantum spin fluctuations, which for frustrated 
systems is a challenging theoretical problem ]|, g|. We 
also show that the low-energy states of our model form 
a large number of valleys separated by energy barriers, 
which prevents the system from reaching thermal equilib- 
rium at low temperatures, but does not result in a spin- 
glass behavior in the absence of quenched disorder j|, || . 

One of the simplest classical frustrated models 
describes the Ising spins <7j = ±1 on sites of a trian- 
gular lattice with nearest-neighbor antiferromagnetic 
interactions: 



E = J^<Ti<Tj, J>0. 



(1) 



This system stays disordered at all nonzero temperatures 
P, ml, while at T = the spins order periodically 
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g. This algebraic order coexists with an extensive zero 
temperature entropy So ~ 0.323fcs per spin q, [j. 

This entropically induced ordering resembles the crys- 
tallization in hard-sphere liquids at high volume fractions 
§, [l0). We show below that the model Eq.Q can indeed 
be mapped on a two-dimensional liquid of topological 
domain walls (kinks), the motion of which is confined to 
one spatial dimension. This mapping also provides useful 
insights into the physics of the quantum version of Eq. (Q) 
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where <r°/' z arc the Pauli matrices and h is the transverse 
field (the classical and quantum frustrated models Eq. (Q) 
and Eq.(H) are referred to below as, respectively, CFM 
and QFM). The critical behavior in the QFM was re- 
cently discussed by Moessner et ah, who argued that the 
ordering of quantum spins is long-ranged at low tempera- 
tures and that the ordered and disordered phases are sep- 
arated by a phase with algebraically decaying spin corre- 
lations || n. Here we show that the difference in the 
critical behavior of the CFM and QFM originates from a 
higher rigidity of the quantum kink crystal. We also per- 
form numerical simulations of the QFM and find an unex- 
pected specific heat anomaly at strong transverse fields. 
Kinks in frustrated Ising model: We shall consider the 
triangular lattice as an array of coupled chains, running 
in the x-direction. As neighboring chains are shifted with 
respect to each other, we will distinguish even and odd 
chains. In each chain we perform the transformation from 
spins to kinks, which are domain walls separating two 
different Neel states and carrying the topological charge 
q = ±1 (kinks and antikinks). The chain energy equals 
2 JN, where N is the number of kinks, independent of 
kink positions. The interchain spin coupling gives rise to 
interactions between pairs of kinks in neighboring chains 
with the potential 



where is the distance between the spins and r] — 1/2 



V(x - x e ) = 2Jg g e sign(x - x e ), 



(4) 



2 





(b) 



FIG. 1: (a) Optimal relative ordering of kinks; the black(whitc) 
circles correspond to kinks with the topological charge +1(— 1); 
(b) Dislocation in the kink crystal with the energy 2J;(c) The kink 
crystal, in which the kinks delocalized over pairs of neighboring 
sites are indicated by dimers; (d) The spin ordering corresponding 
to the kink crystal. 



where q {q e ) and x (x e ) are, respectively, the topological 
charges and ^-coordinates of the kinks in the odd(even) 
chain. This potential only depends on the sign of the 
relative coordinate x a — x e . Therefore, the energy of the 
CFM is completely determined by the number of kinks 
and their relative ordering in neighboring chains. These 
loose interactions between kinks are much different from 
those in unfrustrated models (e.g., the Ising model on 
a square lattice) which grow linearly with the distance 
between kinks, confining them into pairs. 

While in unfrustrated models the density of kinks van- 
ishes at low temperatures, in the CFM kinks are present 
in ground states. In the energetically most favorable 
relative ordering the q = +1 kink in an odd chain has 
the nearest-neighbor q = +1(— 1) kinks in two neighbor- 
ing even chains from the right/left and vice versa (see 
Fig. |l|a), in which case the kink creation energy 2 J is ex- 
actly compensated by the energy of its interactions with 
kinks in neighboring chains, i.e. in the ground states 
kinks cost no energy. For the favorable relative ordering 
the number of kinks N in each chain has to be the same. 
Therefore, the ground states form distinct classes labeled 
by N. Since shifts of kinks that preserve the ordering do 
not change energy, each class still contains a large number 
of states, resulting in an extensive ground state entropy. 

To describe the statistics of kinks in the ground state 
class N, we introduce the 'wave function' ^n{z}, which 
equals the number of the minimal-energy kink configu- 
rations in the lower half-plane for fixed positions of the 
TV kinks in the uppermost chain represented by {z} = 

(zx, Z2, ■ ■ ■ , zn), where zj = e 2m ~ and L x is the chain 
length. One can then add one more chain from above 
and obtain an eigenvalue equation for the wave function 



X N ^ N {z} 



{*'} 



(5) 



where the ~}2' denotes the summation that preserves the 
energetically favorable relative ordering of kinks with 
the coordinates {z} and {z 1 } in two neighboring chains. 
The solution of Eq.(||) is the absolute value of the van 
der Monde determinant 



^n(zi,z 2 , 
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It can also be written in the form of the Slater deter- 
minant of N plane waves e iknXj (n,j — 1,2, ...,N) 
with the wave vectors k n taking values in the Fermi 
sea — kp < k < +kp, where the Fermi wave vector /cf 
is related to the density of kinks n = -P- by kp = nn. 
Thus, Eq.(||) is the ground-state wave function of TV 
fermions in the chain with L x sites and the kinks 
can be identified with the fermions appearing in the 



transfer-matrix solution of the CFM 13 . 

The eigenvalue Xn in Eq.(|5|) is related to the number 
of the ground states in this class, Wn, by Wn = A 
{L y is the number of chains). In the limit L Xl L y 
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and for a fixed kink density n, A at — ^ , ^ n , 

■y L x tan 

where S(n) is the ground-state entropy per site 
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S(n) = - / d<pln ( 2cos n 
7T ./n V 2 
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For n = ^, at which S(n) has its maximum, Eq.(|^) gives 
the value ~ 0.323&B cited above. 

The joint distribution function Pn{{z}) of ./V kinks in 
a chain is obtained by summing over all minimal-energy 
configurations of kinks in the chains both below and 
above this chain, which gives 



P N {z} = *%{z}<x]l\ 



(8) 



The spin correlation function along chains (o~ x o~o) 
= (~ ) x {i") 1 ^^) 7 where K(x) is the number of the 
kinks in the interval [0, at]. Using Eq.(||), the spin 
correlator can be written in the form of the Toeplitz 
determinant: (o- x a ) = (-) x det f nm , where f nrn = 

Snm ~ 7^~ S1 " xf^-m) ■ F° r 1 < s < ii, the determinant 



cx coskpx, so that for kp — ^ we recover Eq.(|]). 
In general, the ground state class with N kinks/chain 
has an algebraic SDW order with the wave vector 
q = 7r(l — n), which we interpret as the 2fcf-instability 
of the kink Fermi sea. For Ising-type models with the 
discrete Zi symmetry such SDW states are very unusual. 
The number of ground states Wn has a sharp peak at 



TV* = L x /3: Wn <x exp 



Though 



-const -j^ (N-N«Y 

the number of classes significantly contributing to the 
CFM partition function stays finite in the thermody- 
namic limit, all of them are essentially identical copies 
of the class with n = 1/3, since the corresponding SDW 



3 



vectors deviate from q = 2tt/3 by an amount oc 1/L X . 
Furthermore, the deviation of the total ground state en- 
tropy per spin S = ln^^ Wjy, where L — L x L y is the 
number of spins, from the entropy of a single class Eq.(0) 
at n — ^ is 0(1/L). Thus, in the thermodynamic limit it 
suffices to consider only one class with 1 kink per 3 sites. 
Quantum model: The transverse field h in Eq.(^) flips 
spins, resulting in hopping of kinks along the chains, 
as well as in creation/annihilation of kink-antikink 
pairs on neighboring chain sites. For h -C J the 
hopping plays the dominant role, as it mixes degener- 
ate classical ground states within each class, whereas 
the kink-antikink pairs cost energy ~ 4J and only 
appear virtually, renormalizing the kink dispersion: 
E(k) ps -2/icosfc + (sin 2 k + \). Here k is the kink 
wave vector, the first term is due to the hopping of 
kinks on neighboring sites, and the second term contains 
contributions from the dressing of kinks and the ground 
state by one virtual kink-antikink pair. For h,T <C J the 
positions of the dressed kinks satisfy the same restrictions 
as in the classical ground states and h is the only relevant 
energy scale. In other words, kinks form a quantum hard 
sphere liquid. The restricted motion of kinks makes the 
quantum system more rigid than the classical one and 
gives rise to phonon-like excitations with velocity oc h. 

In the CFM kinks only crystallize at T = 0. At any 
nonzero temperature the algebraic order is destroyed 
by dislocations in the kink crystal with energy 2J (see 
Fig.|l|b), which in the classical model are unbound. The 
latter leads to the asymptotic behavior of the spin corre- 
lation function at T <C J H , given by the right-hand side 



of Eq.(H) multiplied by e 



where £ = e 2 ^ J is the 



average distance between the dislocations. On the other 
hand, in the QFM phonons result in two-dimensional 
Coulomb interactions U(r) oc hhxr between the dislo- 
cations separated by a distance r, binding them into 
pairs at a finite temperature T% oc h, below which the 
spin correlations decay algebraically |14]. Upon lowering 
temperature, the exponent rj (see Eq.(3)) decreases, as 
the kink crystal becomes more rigid, which ultimately 
leads to the pinning of the kink crystal by the lattice at 
some temperature < T\ , below which the phonons be- 
come gapped and the spin ordering becomes long-ranged 
[ fTo] , |l6| , [ItJ . A cartoon of the quantum kink crystal with 
I kink per 3 sites is shown in Fig.[l]c. To gain kinetic 
energy each kink is delocalized over one bond (such 
bonds are shown as dimers) and the bonds are arranged 
in a way that ensures the energetically favorable ordering 
of kinks. This state corresponds to the SDW state with 
the wave vector q = ^ (see Fig. |I|d), in which the spins 
antifcrromagnetically ordered along the z axis form a 
bipartite hexagonal lattice, while the spins located in the 
centers of the hexagons are oriented along the transverse 
field, since the fields from their neighbors add to zero. 

Thus for weak transverse fields h <C J, the QFM 
describes a kink crystal that melts like a crystal of 
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FIG. 2: The plot of specific heat (circles) and the 'order parameter' 
m (squares) of the QFM vs h/T, for h <C J- The smooth line 
interpolates the part of c(h/T) due to quantum superpositions of 
classical ground states. 



adsorption atoms on a substrate lattice, i.e. the melting 
is preceded by a depinning transition and the solid 
phase is separated from the liquid phase by a 'floating 
crystal phase' with algebraic crystal order |l5|, [l(| p7| . 
In the context of the QFM this phase diagram was 
recently suggested by Moessner et al. Jill who 
studied the critical behavior of the QFM in the vicinity 
the quantum critical point, i.e. in the opposite limit 
of strong transverse fields h ~ J. Two-dimensional 
crystals are also known to melt via a single first-order 
transition due to proliferation of the boundaries between 
degenerate crystal phases 10, Ij], We compare these 
two scenarios to our numerical results. 



Numerical results: We performed Monte Carlo simula- 
tions of the QFM using the continuous time algorithm 
pUf , calculating the temperature dependence of the 
specific heat c and the susceptibilities \q = (Mg)/L, 
where M q = J2j v] cos qxj and q is a multiple of j^- . A 
large value of Xq is a signature of the SDW state with 
the wave vector q. We also calculate the finite system 
'order parameter' m = (M^^) / L. 

The behavior of the model is most spectacular for weak 
transverse fields, h <C J = 1, when one can see a clear 
difference between the classical and quantum regimes. In 
Fig.^ we plot the specific heat and 'order parameter' as a 
function of the ratio T/h, varying over four decades and 
covering both the classical region h T, J and the quan- 
tum region h,T <C J. The 'classical points' are calcu- 
lated at h = 0.01, while the 'quantum data' is a collection 
for many (h, T) points, which fall on smooth curves when 
plotted versus T/h, showing that in the quantum regime 
h is indeed the only relevant energy scale. The specific 
heat has two maxima: one at T ~ J, which also exists in 
the classical model, and another at T ~ h, due to quan- 
tum superpositions of classical ground states. To show 
that the huge degeneracy of the classical model is lifted 
by a transverse field, we calculate the entropy release 
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FIG. 3: The susceptibility x ( a )i the 'order parameter' m (b), 
and the susceptibility x ( c ) versus inverse temperature /3 for h = 
J = 1, = 96 and three different values of L y : 20 (pluses), 40 
(open circles), and 80 (filled circles). Plotted in panel (d) is the 
corresponding temperature dependence of the specific heat. The 
solid line is the specific heat of two-dimensional phonons with a 
gap ~ 0.5. 

related to the low-T maximum, AS — Jq'^-c, where 
h <C T„ <C J- The numerical integration that uses the 
smooth-curve fit of the MC data (see Fig.|J) gives AS = 
0.32fcs, in perfect agreement with the zero temperature 
entropy of the classical model. The 'classical maximum' 
corresponds to the disappearance of defects in the kink 
crystal with the energy ~ J, which in finite systems in- 
duces an algebraic spin ordering with rj = |, just as in the 
CFM at T — 0. The latter is clear from the temperature 
dependence of the 'order parameter', which grows fast at 
T ~ J and stays constant at h < T < J. At T ~ h, the 
'order parameter' grows again, reflecting the increasing 
stiffness of the kink crystal due to the quantum motion of 
kinks, which results in a decrease of the exponent rj and, 
perhaps, in the appearance of the long-range order. The 
smooth temperature dependence of c and m is evidence 
against a first-order transition and is compatible with the 
two phase-transitions scenario, since both transitions are 
expected to be of the Kosterlitz-Thouless type Jl6| . Fur- 
thermore, in a finite system the transition at the upper 
critical temperature T\ ~ h <C J, describing the disloca- 
tion binding does not occur, as all dislocations disappear 
at much higher temperature ~ J. The lower critical tran- 
sition is also rather difficult to identify, since the algebraic 
order with n = 1/9 at T = T 2 + Q is practically in- 
distinguishable from the long-range order at T < T^. 

In Figs. ||(a) and (b) we plot \ = Xitx/z an d rn vs 
temperature, for the strong transverse field h = J = 1, 
L x = 96, and L y = 20,40, and 80. For (3 < 1.5, the sus- 
ceptibility/spin x is independent of L yi corresponding to 
a disordered phase, while for f3 > 2, the 'order parameter' 



m shows little size-dependence for large L y , indicating 
a long-range ordering. On the basis of our data it is 
difficult to conclude whether the intermediate region 
1.5 < /3 < 2 is the 'algebraic' phase with temperature- 
dependent rj or it is a vicinity of a single transition. 
Note, that the specific heat has a rather sharp kink at 
[3 ~ 1.8, at which its behavior changes from approx- 
imately temperature-independent to Independence, 
corresponding to the specific heat of phonons in the 
two-dimensional kink crystal (see Fig. ||(d)). Also the 
susceptibility \ = X^2tt/3 Xq' describing SDWs with 
subdominant harmonics has a peak at (3 — 1.75 (see 
Fig. ||(c)), suggesting a single transition, which may be 
attributed to a sudden loss of rigidity of the kink crystal. 
Topological spin glass?: We found that for (3 J > 2 it is 
effectively impossible to bring a large system into ther- 
mal equilibrium, as it 'freezes' in one of the SDW states 
with q ~ 2ir/3. This 'glassy' behavior is related to the 
existence of different ground state classes, which in the 
quantum case transform into an array of energy valleys 
separated by barriers. At low T the barriers become im- 
penetrable, as the tunneling between neighboring classes 
requires creating/annihilating a kink-antikink pair in all 
chains. Does such an energy landscape lead to a spin- 
glass behavior in the absence of quenched disorder and 
should the 'kink' in c(T) at T - 0.6 (see Fig. |c) be in- 
terpreted as a spin-glass transition? We believe that the 
answer is negative, because (as in the CFM) the number 
of important energy minima in the thermodynamics limit 
stays finite, all of them describing essentially the same 
state. In our low-T simulations we select the state with 
q = 2ir/3 by an appropriate initial spin configuration. 

In conclusion, we explained unusual properties of the 
classical and quantum frustrated Ising models by map- 
ping them on a system of kinks, which behave like hard 
spheres. We showed that the quantum spin dynamics 
makes the kink crystal more rigid, resulting in a long 
range ordering at low temperatures. We studied numer- 
ically the temperature dependence of specific heat and 
spin susceptibility in weak and strong transverse fields. 
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